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We report on the growth and electrical characterization of a series of two-dimensional hole systems 
(2DHSs) used to study the density dependence of low temperature mobility in 20 nm GaAs/AlGaAs 
quantum wells. The hole density was controlled by changing the Al mole fraction and the setback 
of the delta-doping layer. We varied the density over a range from 1.8 x 10^° cm~^ to 1.9 x lO'^^ 
cm~ finding a nonmonotonic dependence of mobility on density at T = 0.3 K. Surprisingly, a peak 
mobility of 2.3 x lO'' cm'^/Vs was measured at a density of 6.5 x 10^'^ cm^'^ with further increase 
in density resulting in reduced mobility. We discuss possible mechanisms leading to the observed 
non-monotonic density dependence of the mobility. Relying solely on interface roughness scattering 
to explain the observed drop in mobility at high density requires roughness parameters which are 
not consistent with measurements of similar electron structures. This leaves open the possibility of 
contributions from other scattering mechanisms at high density. 



Two-dimensional hole systems (2DHSs) on (001) ori- 
ented GaAs offer an interesting alternative to the 
more widely studied two-dimensional electron systems 
(2DESs). 2DHSs on (001) GaAs liave effective masses 
roughly 4.5 to 7.5 times larger—"— than that in cor- 
responding 2DESs which increases the importance of 
Coulomb interactions relative to the kinetic energy re- 
sulting in enhancement of importance of many-body ef- 
fects. In addition, the p-wave symmetry of the valence 
band in GaAs leads to a much reduced hyperfine coupling 
of hole spins to the atomic nuclei which makes them an 
exciting alternative to electrons for quantum dot spin- 
based qubits^"— . The presence of spin-orbit coupling and 
light/heavy hole mixing in the valence band of GaAs also 
allows extensive band structure engineering^"— . This fea- 
ture has been exploited to alter the nature of ground- 
states in the quantum Hall regimei ^°d^ 

Here we describe our efforts to understand the limits 
to low temperature mobility for (001) 2DHSs. Contin- 
ued improvement in 2DHS quality is motivated by the 
well-established paradigm for 2DESs that increased low- 
temperature mobility often leads to the observation of 
new correlated groundstates^^. Historically, improve- 
ment to the low temperature mobility of 2DHSs has 
lagged behind that of 2DESs due to the lack of a p- 
type dopant in GaAs that docs not diffuse or segregate 
significantly at typical molecular beam epitaxy (MBE) 
growth temperatures ~ 635 °C. Si can act as a low- 
diffusivity acceptor on the (311)A face of GaAs, but sub- 
sequent transport experiments are known to be compli- 
cated by a significant mobility anisotropy due to sur- 
face corrugatioriH. However, recent use of low diffusiv- 
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ity carbon doping (C-doping)^"— has rapidly led to low 
temperature mobilities > 10^ cm^/Vs without the ac- 
companying transport anisotropy. Purely from a growth 
standpoint, then, there does not appear to be any reason 
why low temperature hole mobilities should not approach 
that of electrons once scaled by the appropriate effective 
mass. Presently it is widely believed that uniformly dis- 
tributed ionized background impurities limit the mobility 
in the best 2DES9i^. However, the highest hole mobility 
reported to datoi^ of 2.6 x 10^ cm^/Vs is still about a 
factor of two lower than record mobility 2DESs grown in 
the same MBE chamber— once the heavy hole to electron 
effective mass ratio of 0.4me : 0.067me is taken into ac- 
count. The question then remains, if sufficiently reducing 
background impurities^^ is the main obstacle to reaching 
an electron mobility of 100 x 10® cm^/Vs, what are the 
key ingredients to a hole mobility of 15 x 10® cm^/Vs? 

In order to answer this question, we have begun to 
explore the impact of varying structural parameters on 
the resulting mobility. Samples in this work were grown 
in a customized Veeco Genii MBE which has recently 
achieved electron mobilities > 20 x 10® cm^/Vs and ex- 
tremely large excitation gaps for the fragile v=b/2 frac- 
tional quantum Hall state. C-doping was performed with 
a carbon filament capable of producing a doping rate of 
2.8 X IQ^^ cm~^/sec at a total power (including para- 
sitic dissipation) of ~ 150 WJ^ In this experiment, we 
utilized a 20 nm quantum well situated 190 nm below the 
surface and asymmetrically i5-doped from above at a set- 
back d of 80; 110; or 140 nm. The Al mole fraction x was 
also varied between 0.07 and 0.45 to allow further tuning 
of the 2DHS density. Table U summarizes the structures 
grown in the experiment, and Fig. [1] shows the epilayer 
design. Square samples were prepared using InZn con- 
tacts annealed at 430 °C for 15 minutes in H2/N2 forming 
gas. Characterization was performed in the dark and af- 



TABLE I: Summary of structural parameters including S- 
doping setback distance d, Al mole fraction around the 
dopants Xd, Al mole fraction surrounding the quantum well 
Xji,, 2DHS density p, and T = 300 mK mobility after illumi- 
nation fi. 
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FIG. 1: Layer structure of devices in this experiment. Note 
the use of two different Al mole fractions Xw and Xd in some 
of the devices as indicated in Table HI 



ter illumination with a red LED at T = 300 mK using 
standard lock-in techniques, and the density was deter- 
mined from quantum Hall effect (QHE) minima. Illumi- 
nation typically resulted in ~ 3-5% increase in density 
and as much as a 27% increase in mobility for low den- 
sity samples. Transport data also showed a qualitative 
improvement after illumination, indicating that illumi- 
nation increases the homogeneity of the 2DHS and has 
a favorable impact on the screened disorder potential. 
Figure [5] shows transport data of the highest mobility 
sample and a low density sample; the number of nascent 
fractional QHE features attest to the sample quality. 
Fig. [3^ shows the measured mobility as a function of 
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FIG. 2: (Color online) Magnetotransport a.t T = 300 mK after 
illumination with a red LED of (a) peak mobility sample and 
(b) low density sample that exhibits many nascent fractional 
QHE features. 



density for various values of d. We note that remote ion- 
ized impurity (RI) scattering does not play a significant 
factor in limiting the mobility since within experimen- 
tal uncertainty there is no meaningful difference between 
the mobility at different values of d for the same den- 
sity. However, increased d should allow these samples 
to be gated to ultra-low densities before RI scattering 
begins to cause the mobility to rapidly drop off with fur- 
ther decreased densitjJ^. The most interesting feature 
of the data in Fig. [3^, however, is the strongly non- 
monotonic dependence of the mobility on density. For 
2DESs in this density range with such a large value of d, 
one would expect the mobility to monotonically increase 
with densitj*i^'iii^'2i following a power law dependence 
fi (X p°' where a '^ 0.6 - 0.8 with ionized background 
impurity (BI) scattering being the dominant scattering 
mechanism. In analyzing our results, we first note that 
the eifective mass is known to vary throughout the den- 
sity range of our samples due to the valence band non- 
parabolicity arising from light- and heavy-hole band mix- 
ing. By performing a linear fit to cyclotron resonance 
data on 2DHSs in (001) 20 nm quantum wells in refs. 
[lIQ and assuming the cyclotron mass plateaus at 0.5me 
at high density we estimate the transport lifetime for our 
structures as shown in Fig. [SJd. The transport lifetime, 
however, follows the same non-monotonic behavior as the 
mobility which indicates a competition between different 
scattering mechanisms throughout the density range of 
our experiment in addition to the changing mass. 
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FIG. 3: (Color online) (a) T = 300 mK mobility after illu- 
mination with a red LED as a function of density for various 
dopant setback distances d. Solid lines are guides to the eye. 
For fixed d the density was controlled by varying the Al mole 
fraction x. Samples were grown in random order to avoid con- 
tinued machine clean-up from skewing the observed trend in 
mobility. Samples A-D were grown with varying x at fixed p 
to test the effect of alloy and interface roughness scattering on 
the mobility (see text), (b) Transport lifetime estimated as a 
function of density. Inset: Effective mass for our structures 
as a function of density extrapolated from refs. [l|[ 



To shed further light on possible scattering mecha- 
nisms, we have performed a series of scattering calcula- 
tions including the effects of BI, RI, alloy, and interface 
roughness (IR) scattering. We follow the derivation of the 
transport relaxation time in [23 which assumes T" = 
and neglects intersubband scattering, multiple scattering 
events, and correlation between ionized impurities. This 
simple calculation is intended to elucidate the expected 
trend of the mobility as the density is increased and de- 
termine if scattering mechanisms dominant in 2DESs can 
qualitatively explain our observations. More sophisti- 
cated calculations have been made by S. Das Sarma and 
coworkers ii2i2^i2i Transport relaxation times are calcu- 
lated individually and then the total mobility is calcu- 
lated using Mathicsscn's rule. For BI and RI scattering 
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where m* is the hole effective mass (as estimated in Fig. 
|3b), h is the reduced Planck constant, is the scatter- 
ing angle, ZiCi is the impurity charge, e is the dielectric 
constant of the semiconductor, q is the scattering vector, 
kp is the Fermi wavevector, qxr is the Thomas-Fermi 
screening wavevector, Ni{z) is the i*^ impurity distribu- 
tion, and the form factors are given by 

5.(q) = / x\z)x\z')eM-<l\z - ;^'|)dzdz' (2) 



(3) 



where x(z) is the self-consistcntly calculated^ envelope 
function in the effective mass approximation. For the 
BI calculation we use a three-dimensional impurity con- 
centration N'^D as a fitting parameter and find the best 
agreement with the experimental data for N^£> = 2 x 
10^'^ cm""^. We use a remote impurity sheet concentra- 
tion Nri equal to the hole concentration p. A more real- 
istic value of Nri could also include some of the ionized 
impurities due to the surface compensation; however, we 
assume a simple parallel-plate capacitor model of the 
surface-(5-layer charge and thus neglect the surface com- 
pensation contribution to Nj^j. This neglect of charge 
due to surface compensation is typical in these types of 
calculations fi2i2£i2i For our purposes, though, the exact 
value of Nui is not important since it will not change the 
qualitative dependence of the RI- limited mobility as p is 
varied. 

To calculate alloy scattering we use the virtual crystal 
approximation with a square well potential limited over 
a spherical range^^ which is independent of temperature 
in 2D systems^S. The alloy limited relaxation lifetime 
is unscreened due to its short range nature and given 
by^i^ 
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where a = 0.565 nm is the lattice constant of the com- 
pound semiconductor, 17 is the volume of the scattering 
potential given by il = (4/3)7rr'^, and r = (\/3/4)a is the 
nearest-neighbor separation. There is a broad range of 
estimates of the magnitude of the scattering potential U 
in the literature^, ranging from 0.12 to 1.56 eV. We take 
U~ 1 eV (as suggested in ref. [23l) as a rough estimate. 
To examine the possible effect of interface roughness 
scattering, we employ a simple model which makes use 
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FIG. 4: (Color online) Dashed lines show a comparison of the 
self-consistently calculated valence band edges (dashed lines) 
for the high density sample # 3 and and a single heterojunc- 
tion sample with x = 0.45. Solid lines show a comparison of 
the self-consistently calculated wavefunction for sample # 3 
and the Fang-Howard variational wavefunction. 



of the Fang- Howard variational wavefunction and associ- 
ated potential^ which takes the distortion of the wave- 
function with increased density into account. In this 
model, the IR scattering rate is given by^>^ 
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where the wavefunction used to calculate g^ is 
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We take one monolayer roughness height to be a reason- 
able estimate and thus set A = 0.2825 nm and use A as 
a fitting parameter with the result that A = 6 nm. 

As a justification for using the Fang- Howard wavefunc- 
tion to model our asymmetric quantum well system, we 
show in Fig. 2] a comparison of the self-consistently cal- 
culated valence band edges for quantum well (QW) sam- 
ple # 3 and a single heterojunction (SHJ) along with 
the self-consistently calculated wavefunction for the QW 
structure and the Fang-Howard wavefunction which is 
often taken as an approximation of the wavefunction in 
SHJ structures. The band edges show that the bottom 
barrier of the QW changes the confining potential very 
little, and the high density samples (where IR scattering 
could be important) can therefore be approximated by 
the Fang-Howard model. 

The results of our calculations are compared with the 
d = 80 nm experimental results in Fig. \5\ It is clear that 



even with a changing effective mass and wave function 
profile the BI and RI limited mobilities steadily increase 
with increasing density and therefore cannot account for 
the drop in mobility at high density. The exact contri- 
butions of alloy and interface roughness scattering are 
initially less clear, however. We will address alloy disor- 
der first. 

If U is large enough, alloy scattering could conceivably 
contribute to the drop in mobility seen in the experimen- 
tal data. Before continuing, it should be noted that the 
slight increase in the calculated alloy-limited mobility at 
high density is simply due to the saturation of the ef- 
fective mass as shown in the inset of Fig. |3}d. We have 
repeated the calculations (not shown) without forcing the 
mass to plateau at 0.5me, but even with a mass as high 
as O.TrUe at high density the alloy limited mobility does 
not appear to be limiting the total mobility. To test the 
contribution of alloy scattering we grew a series of four 
test structures (labeled A-D in Fig. [3^) in which Xd, 
the Al concentration starting 25 nm above the quantum 
well (e.g. around the (5-doping layer), was kept fixed to 
keep the density constant while x^j, the Al concentra- 
tion around the quantum well, was varied between x^ = 
0.16 and x^, = 0.24. Samples A, B, and 3 (xd = 0.45) 
suggest that Xw has no impact on the mobility, though 
there is scatter in the resulting density which we attribute 
to wafer-to-wafer variation and possible variation in the 
illumination. Samples C, D, and 5 (xd = 0.35), how- 
ever, suggest that increased x^, does cause the mobility 
to decrease somewhat. Most importantly, this is the op- 
posite trend that would be expected if alloy scattering 
per se were limiting the mobility. Our calculations for 
the test structures (not shown) and ref. [33| predict that 
the alloy-limited mobility would increase for increased Xw 
since as x^^j is increased for fixed density the wavefunc- 
tion is more confined. This in turn causes the integral 
of x"* to decrease faster than the x{l — x) term increases 
in equation ((4|) resulting in a decrease in the alloy scat- 
tering rate for increased Xw ■ The results from this set of 
structures is consistent, however, with the theory that Al 
getters impuritiesi^, thus an increase in Xw would locally 
increase Nbi and the associated scattering. Regardless, 
samples C, D, and 5 show that the negative side effects of 
increasing x^j arc not enough to explain the data of Fig. 
[3l If the increase in x^ was dominating the mobility 
we would expect test structures A and C to have signif- 
icantly higher mobilities than the peak mobility sample 
#11 due to the higher hole density of the test structures. 
This, however, is clearly not the case. 

Finally, our fit seems initially to indicate that IR scat- 
tering is limiting the mobility at high density. However, 
whenever parameters can be freely adjusted caution must 
of course be exercised to obtain physically meaningful re- 
sults. The dashed pink line in Fig. [5] shows the IR lim- 
ited mobility for a SHJ 2DES in the Fang-Howard model 
while the pink star shows a 2DES SHJ structure with x 
= 0.35 grown during the course of this experiment. Ev- 
idently the Fang-Howard calculation overestimates the 
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FIG. 5: (Color online) Comparison of d = 80 nm experimental 
data with mobility calculations. Nbi and A are used as free 
parameters to obtain a good fit to data. Pink star represents 
SHJ 2DES grown during this experiment. 
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FIG. 6: (Color online) Comparison of d = 80 nm experimental 
data with mobility calculations. A and A are varied in order 
for mobility calculations to obtain agreement with an electron 
structure grown during this experiment (see text). 



IR scattering by at least a factor of four. Repeating our 
self-consistent calculation for Bl, Rl, and alloy scattering 
in this 2DES SHJ using the impurity concentrations and 
alloy potential listed in the inset of Fig. [5] we find that 
the IR-limited mobility at a density of 2.4 x 10^^ cm~^ 
would have to be 86 x 10^ cm^/Vs to fit the measured 
total mobility of 7.9 x 10^ cm^/Vs. To get such a high 
IR- limited mobility we are forced to set A = 0.1 nm and 
A = 2.2 nm. Figure [5] shows the result of our calcula- 
tion for the hole structures using these smaller roughness 
parameters. With these reduced roughness parameters 
there is no longer a good fit to the hole data at high den- 
sity as the IR term makes almost no contribution to the 
total mobility, though we still obtain a good fit at low 
to medium density. We conclude that our crude model 
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FIG. 7: (Color online) Shubnikov-de Haas oscillations of high 
density d = 80 nm, x = 0.45 device. Inset (a): Sketch of 
the spin-split heavy hole and light hole ground states in a 
quantum well. Inset (b): Index of extrema in Kxx vs. B~^. 



The high field slope gives the total density of 1.8 x 10 cm~ , 
and the low field slope gives the lighter sub-band density of 7 
x 10^" cm~^ while the difference in the two gives the second 
sub-band density of 1.1 x 10^^ cm~^. 



of interface roughness scattering cannot simultaneously 
account for our experimental data in both electrons and 
holes and are thus hesitant to conclude that interface 
roughness scattering is the dominant source of our drop 
in mobility at high density. Similar discrepancies between 
electron and hole data have been noted in ref. j35| . 

Another possible scattering mechanism that must be 
kept in mind at high density is scattering between the 
electric subbands of the quantum well which is known 
to degrade the mobility in high density 2DESsi^ To es- 
timate the possibility of such scattering, we use a finite 
square well with a barrier height of 230 meV and an ef- 
fective mass of 0.5me, which corresponds to our highest 
density sample. This estimate results in an energy spac- 
ing of 5.0 meV between the heavy hole ground and first 
excited state. Assuming a light hole mass along the (001) 
directional^ of 0.090me the spacing between the heavy 
hole and light hole ground states is 6.4 meV. In both 
cases this energy spacing is significantly larger than the 
Fermi energy Ep — ^ ~ 0.9 meV which precludes 
a significant contribution from intersubband scattering 
between the electric subbands. 

Next, we note the presence of beating in the 
Shubnikov-de Haas oscillations in Fig. [7] which is indica- 
tive of B = spin-splitting. Such spin-splitting is known 
to occur in structurally-asymmetric devices^'^'^ due to 
Rashba spin-orbit coupling^. Wc sketch the qualitative 
effect of this splitting in inset (a) of Fig. [T] As the 2DHS 
density is increased, the electric field (and hence spin 
splitting) in the well is also increased. Furthermore, it is 
know that the presence of a parallel channel can result 
in a Hall density different from the sum of the subband 
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FIG. 8; (Color online) (a) Estimate of the Hall density ex- 
pected from Eq. [7] using the measured subband densities 
pi = 1.1 X 10" cm"^ and p2 = 7 x 10^° cm"^. (b) Estimate of 
the expected measured mobility \iHaU if the high density sub- 
band pi — 1.1 X 10^^ cm^^ has a high mobility pii = 2 x 10® 
cm^/Vs. For the second subband density (dashed red line) 
P2 = 7 X 10^" cm~^ measured in Fig. [7] the expected mea- 
sured mobility i^Haii > 1-75 x 10® cm^/Vs. 



densities and a measured mobility different from that 
of either subband even in the absence of intersubband 
scattering. In our case, wc assume that the two parallel 
channels are non-interacting B = spin-split subbands of 
the heavy hole ground state. The measured Hall density 
Vtiaii and mobility ^.uaii in the absence of intersubband 
scattering are given by^ 
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where y^i(2) and Pi{2) arc the mobility and density, re- 
spectively, of the first (second) subband. Figure [5^ il- 
lustrates the Hall density as a function of the subband 
mobilities in our peak density sample (sample # 3) pre- 
dicted by Eq. [7] using the subband densities extracted in 
Fig. [71 It is clear from Fig. [5^ that in order to measure 
a Hall density ~ 1.8 x 10^^ cm~^ the subband mobilities 
should be comparable, though the high density subband 



should have a slightly higher mobility. In order to esti- 
mate the effect of the presence of two subbands on the 
measured mobility, we therefore assume that the high 
density subband is also the high mobility subband. In 
order to determine if the presence of the lower mobility 
subband could by itself account for the drop in mobility 
seen in Fig. |31 we assume that the high mobility sub- 
band is unchanged from the peak total mobility value 
(^ 2 x 10^ cm^/Vs) at low density. Figure ^p shows 
what we would thus expect to measure as a function of 
density and mobility in the low mobility subband if the 
high mobility subband has a density pi = 1.1 x 10"'^^ cm~^ 
as we measure in Fig. [T] For the measured second sub- 
band density of p2 = 7 x 10^° cm^^ (dashed red line) we 
see that this parallel subband effect would not decrease 
the measured mobility below ^ 1.75 x 10^ cm^/Vs. We 
therefore conclude that the presence of a second, possibly 
low mobility B = spin-split subband cannot explain our 
observed drop in mobility at high density in the absence 
of intersubband scattering. 

A final possible mechanism for the observed drop in 
mobility at high density is intersubband scattering be- 
tween the spin-split subbands of the heavy hole ground 
state of the quantum well. The question remains, how- 
ever, whether or not there exists a potential capable 
of coupling the spin-split sub-bands and causing back- 
scattering. Such scattering is typically neglected in the- 
oretical calculations of the mobility due to the assumed 
lack of a significant spin-fiip mechanism^, though in- 
tersubband hole-hole scattering in inversion-asymmetric 
structures is not without precedenti^ At this time more 
theoretical work is needed to resolve the relative contri- 
butions of the different scattering mechanisms. 

In conclusion, we have performed an experimental 
study of the density dependence of mobility in C- 
doped (001) GaAs/AlGaAs quantum wells by varying the 
dopant setback d and Al mole fraction x. The mobility 
was seen to depend non-monotonically on the density. At 
low density the mobility increased with density. The T 
= 300 mK mobility was found to peak at a value of 2.3 x 
10^ cmVVs at a density of 6.5 x 10^° cm-^. This 2DHS 
mobility is among the highest ever reported. Increasing 
the density further, however, resulted in a sharp drop 
in mobility. Scattering calculations indicate that back- 
ground ionized impurities and remote ionized impurities 
will not lead to a decrease in mobility at high density 
even with a changing effective mass, and alloy scattering 
cannot account for all of our experimental results from 
various test structures. Interface roughness scattering 
contributions remain unclear due to the difficulty in ob- 
taining physically reasonable roughness parameters that 
predict both electron and hole mobilities. Beating in the 
Shubnikov-de Haas oscillations in our high density sam- 
ples is indicative of zero-field spin-splitting which leaves 
open the possibility of an intersubband scattering con- 
tribution to the mobility. Further theoretical work is 
needed to determine the mechanism and magnitude of 
such a contribution. 
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